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Abstract 

There are several lines of evidence indicating that the ultra-relativistic out- 
flows powering gamma-ray bursts (GRBs) are collimated into narrow jets. 
However, these are indirect, and the jet structure is rather poorly con- 
strained. What is more, the jet dynamics have still not been investigated 
in detail. It has been suggested that the observed variety between different 
long duration events, ranging from bright spectrally hard GRBs, to dimmer 
and spectrally softer X-ray flashes (XRFs) may be largely due to different 
viewing angles (or lines of sight) relative to rather similar relativistic jets. 
Here we describe the current state of knowledge on these topics, explain 
some of the most relevant physics behind some of the basic principles, and 
discuss prospects for the future. 

11.1 Evidence for bulk relativistic motion in gamma-ray bursts 

The first line of evidence for ultra-relativistic bulk motion of the outflows 
that produce GRBs arises from the compactness argument. It relies on the 
observed short and intense pulses of gamma rays and their non-thermal en- 
ergy spectrum that often extends up to high photon energies. Together, 
these facts imply that the emitting region must be moving relativistically. 
In order to understand this better, let us first consider a source that is either 
at rest or moves at a Newtonian velocity, /3 = u/c <C 1, corresponding to a 
bulk Lorentz factor F = (1 — /3^)~^/^ ~ 1. For such a source the observed 
variability timescale (e.g., the width of the observed pulses) At, implies a 
typical source size or radius R < cAt, due to light time travel effects (for sim- 
plicity we ignore here cosmological effects, such as redshift or time dilation). 
GRBs often show significant variability down to millisecond timescales, im- 
plying R < 3 X 10'^(At/ 1ms) cm. At cosmological distances their isotropic 
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equivalent luminosity, L, is typically in the range of lO^'^ — lO^^ergs"^. 
In addition, the (observed part of the) eF^ GRB spectrum typically peaks 
around a dimensionless photon energy of e = -Eph/mgC^ ~ 1, so that (for 
a Newtonian source) a good fraction of the total radiated energy is carried 
by photons that can pair produce with other photons of similar energy. [F 
is the radiative flux and -F^ = dF/de). A simple estimate of the opacity to 
pair production (77 — )• e^e~) usually results in a huge optical depth for this 
process, T^-y{e) ~ a'Tn^{l/ e)R ~ a^Lij^/A'K'meC'R^cr'TLiii./A'KmeC^^t ~ 
10i^(Li/J105iergs-^)(lms/At), where = dL/de and (Tt is the Thom- 
son cross section (Granot et al. 2008). Such huge optical depths are clearly 
inconsistent with the non-thermal GRB spectrum, which has a significant 
power-law high-energy tail. This is known as the compactness problem (Ru- 
derman 1975). 

If the source is moving relativistically toward us with a bulk Lorentz factor 
r ^ 1, then in its own rest frame (where quantities measured in that rest 
frame are denoted by a prime) the photons have much smaller energies, e' ~ 
e/r, while in the lab frame (i.e., the rest frame of the central source) most of 
the photons propagate at angles ^ l/F relative to its direction of motion (see 
Sect. lll.2] on aberration of light). The latter implies that in the lab frame the 
typical angle between the directions of the interacting photons is 612 ~ 1/r, 
which has the following consequences. First, it increases the threshold for 
pair production, 8162 > 2/(1 — cos ^12), to £162 ~ (compared to e'ie'2 ^ 1 
for the roughly isotropic distribution of angles between the directions of the 
interacting photons in the rest frame of the source, where 9'i2 ~ 1)- Thus, 
Li/e needs to be replaced by Lr^/e = r^(^""^Li/g, where LqC^"" at 

high photon energies (corresponding to dNpi^/de oc i.e., a is the high- 
energy photon index). This results in an additional factor of r^(^~") in the 
expression for Tjj{e). Second, the expression for the optical depth includes a 
factor of 1— cos 612 (that represents the rate at which photons pass each other 
and have an opportunity to interact) which for a stationary source is ~ 1, 
but for a relativistic source moving toward us is ~ T~^. Finally, the emission 
radius can be as large as R^F^cAt, which introduces an additional factor 
of ~ in the expression for r^^. Altogether, T-y-y(e) is reduced by a factor 
of ~ r2(°+i) (i.e., cx r-2("+i)). Since typically a ~ 2 - 3, this usually 
implies F > 10^ in order to have r^^ < 1 and overcome the compactness 
problem. Using similar arguments, the lack of a high-energy cutoff due to 
pair production in the observed spectrum of the prompt 7-ray emission in 
GRBs has been used to place lower limits on the Lorentz factor of the outflow 
(Krolik k Pier 1991, Fenimore et al. 1993, Woods & Loeb 1995, Baring & 
Harding 1997, Lithwick &: Sari 2001). Recently, high-energy observations 
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Fig. 11.1. The temporal evolution of the radio afterglow image size of GRB 030329 
{left panel) and the implied average apparent expansion velocity in units of c {right 
panel). (Taken from Pihlstrom et al. 2007). 



by the Fermi Gamma-Ray Space Telescope have enabled particularly strict 
(r> 10^) limits to be placed for both long (Abdo et al. 2009a,b) and short 
(Ackermann et al. 2010) GRBs. 

A different, and somewhat complementary, line of evidence for rclativistic 
bulk motion in GRBs comes from estimates of the afterglow image size at 
relatively late times (weeks to years) in the radio. The afterglow image size 
increases with time, and therefore it can be (marginally) angularly resolved 
only at late times, and only for relatively nearby and radio bright afterglows. 
At such late times the flow is much less rclativistic, with a much more 
modest Lorentz factor. The size of the afterglow image at a single epoch 
can be estimated from the quenching of diffractive scintillation in the radio 
afterglow (Goodman 1997), as the angular size of the source becomes larger 
than that of the relevant density fluctuations in the interstellar medium of 
our Galaxy, which is estimated to be roughly ~ 3 nas. When the source 
angular size Og is smaller than Od then ~ (^d/^s)^ ^ Ij and different sub- 
images are formed, which produce a random diffraction pattern resulting in 
frequency dependent brightness fluctuations of order unity {AF/F 1). For 
GRB 970508 this implied an apparent source size of ~ 10^^ cm at t ~ 30 days, 
or an average apparent expansion velocity close to the speed of light, c, 
during the first month (Frail et al. 1997, Taylor et al. 1997, Waxman et al. 
1998). The flux below the synchrotron self-absorption frequency can also be 
used to constrain the size of the emitting region (e.g., Katz & Piran 1997, 
Granot et al. 2005). 

A more direct measurement of the image size, as well as its temporal evolu- 
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tion, can be obtained through very long base-hne interferometric techniques 
at radio wavelengths using, e.g., the Very Long Baseline Array (VLB A). 
This has been possible for only one radio afterglow so far - GRB 030329 at 
a redshift ol z = 0.1685 (Taylor et al. 2004, 2005, Pihlstrom et al. 2007), 
since it requires a bright, relatively nearby event (2:^0.2). Nevertheless, 
the measured source size and its temporal evolution imply a relativistic and 
decelerating apparent expansion velocity (Fig. Ill.ip . with an average value 
of ~ 6c over the first 25 days, and transition to a sub-relativistic expansion 
after about one year. 



11.2 Aberration of light and the Doppler effect 

In relativistic sources such as GRBs, aberration of light (also known as 
relativistic beaming) and the Doppler effect play an important role. These 
two effects refer to the change in direction and frequency, respectively, of 
electromagnetic waves (or particles - photons) between different frames of 
reference - in our case between the rest frame of the emitting fluid (that 
is referred to as the comoving frame), and the lab frame (the rest frame of 
the central source in which the external medium is at rest and the emitting 
jet material is moving relativistically) . These effects can be easily derived 
from the Lorentz transformation of the 4- vectors of an electromagnetic wave 
k^^ = [uj/c, k) or from the photon energy-momentum P^^ = fik^ = E{l/c, k). 
Some more intuition may be gained by the following derivations using the 
Lorentz transformation: t = T{t' + f3x' /c), x = T{x' + f3ct'), y = y', z = z' . 
This implies 

_ dx dx' + Pcdt' v'^ + Pc 



dt dt' + pdx'/c l + ^v'^/c 



(11.1) 



For a photon, Vx/c = cos 6* = fi and v'^/c = cos9' = fi' are the cosines of 
the angle between its direction of motion and the x-direction (defined as the 
direction of motion of the primed frame relative to the un-primed frame) in 
the two rest frames, respectively. Therefore, 

where the second equality is derived from symmetry (interchanging primed 
and un-primed quantities and replacing /3 by — /3). This is the formula for 
aberration of light (or relativistic beaming). It is demonstrated in Fig. 111.21 
for a point source emitting isotropically in its own rest frame. For a rela- 
tivistic source moving with a Lorentz factor T = (1 — /3^)~^/^ ^ 1 (in the 
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Fig. 11.2. Aberration of light: the arrows show the directions of photons in the lab 
frame for a point source that emits isotropically in its own rest frame and moves to 
the right at different values of the four- velocity r/3. For F ^ 1, half of the photons 
(and 3/4 of the radiated energy) are within an angle of 1/F around the source's 
direction of motion (between the grey arrows, which correspond to 0' — 90°). 

lab frame), half of the photons and most (3/4) of the emitted energy are 
within an angle of l/F around its direction of motion. 

The formula for the Doppler factor, T) = v/v' , can be derived noticing that 
the phase of an electromagnetic wave must be Lorentz invariant, since if the 
electric and magnetic fields vanish in one frame then they must vanish in all 
frames. In particular, over one period of oscillation of the electromagnetic 
field, P = l/i^ = A/c, a photon travels a distance equal to its wavelength 
A, and similarly P' = = A'/c, so that V = A'/A. Note that because 
of this simple consideration A must be equal to the difference in the path 
length to the observer (measured simultaneously in the lab frame) between 
two hypothetical photons emitted with a single period time difference (i.e., 
c times the difference in the arrival time to the observer of such hypothetical 
photons). Over a single period, a time P' = A'/c elapses in the comoving 
frame, which due to time dilation corresponds to a time Atiab = TA'/c in the 
lab frame over which the source travels a distance of I = /3cAtiab = r/3 A'. 
During the same time the first hypothetical photon has traveled a distance 
of cAtiab = TA', which implies, from simple geometrical arguments, that 
\ = T\' -1^1 = A'r(l - /3/u) and therefore V = A'/A = l/r(l - |3^l). 

Due to symmetry considerations similar to those mentioned above, V = 
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u' /u = l/r(l + /3/i') and therefore V = l/V = T{1 + In particular, 

l/r(l — j3^) = r(l + /3^'), which is equivalent to eq. 111.21 Expressing T> 
in the comoving frame makes it easier to calculate various quantities for 
an isotropic emitter in this frame, which is broadly expected in most cases 
of astrophysical interest. For example, it is easy to see that the average 
Doppler factor for such an isotropic emitter is simply F, while a fraction 
(2 + /3)/4 of the energy in the lab frame is carried by photons with 9' < 90° 
that correspond to cos 6* > /3 (or < 1/F for F ^> 1). 

Finally, we note that due to symmetry (no preferred direction in the plane 
normal to /3 = v/c), the azimuthal angle of the photon direction of motion 
is equal in the two frames, (j) = (j)', and thus the differential solid angles 
dQ = d(j)dfj, and dO,' = d4>'dii' are related by dQ' /dO, = dfi' /dfj, = V^, using 
ea. fTOl 

11.3 Evidence for jets in gamma-ray bursts 

In contrast to Active Galactic Nuclei (AGN) and Galactic micro-quasars, in 
many of which the relativistic jets are well resolved and their structure can 
be studied directly, the images of GRBs are usually unresolved, or at best 
the image of their late time radio afterglow is only marginally resolved (i.e., 
the image size is still somewhat smaller than the instrumental beam size, 
and its size and shape can be only rather crudely estimated). Therefore, 
the evidence for jets (i.e., highly collimated outflows) in GRBs is mainly 
indirect. The main arguments, or lines of evidence, in favor of jets in GRBs 
are as follows. 

First, other known sources of relativistic outflows (such as micro-quasars 
and AGN) are collimated into narrow jets, so it is natural to expect a sim- 
ilar behavior for the relativistic outflow in GRBs, if indeed the underlying 
process for launching jets are similar (e.g., accretion onto a black hole). 

Second, the very high values deduced for the energy output in gamma rays 
assuming isotropic emission, are inferred for GRBs with known 

redshifts, approach and in some cases even exceed Mqc^ (for GRB 080916C 
-£'7,iso ~ 4.9M0C^; Abdo et al. 2009a). Such extreme energies in an ultra- 
relativistic outflow are hard to envision in models involving stellar mass 
progenitors. If the outflow is collimated into a narrow jet (or bipolar jets) 
that occupies a small fraction <C 1 of the total solid angle, then the strong 
relativistic beaming due to the very high initial Lorentz factor (Fq ^ 100) 
will cause the emitted gamma rays to be similarly collimated. This reduces 
the true energy output in gamma rays by a factor of f^^ to = fbE-y^iso, 
thus significantly reducing the energy requirements (Rhoads 1997, Halpern 
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et al. 1999, Sari et al. 1999). In addition, there is good (spectroscopic) 
evidence that at least some GRBs of the long-soft class (Kouveliotou et al. 
1993) are coincident (to within a few days) with a core collapse supernova 
belonging to the Type Ic category (Stanek et al. 2003, Hjorth et al. 2003). 
In such cases the average Lorentz factor must be (F) < 2 for a spherical 
explosion, since the accreted mass is not expected to significantly exceed 
the ejected mass, and only a fraction of the rest energy of the former can 
provide the kinetic energy for the latter. Therefore, only a minute fraction 
of the ejected mass can reach T > 100 that is required in order to power the 
GRB. Hydrodynamic analysis of spherical blastwaves (Tan et al. 2001, Perna 
& Vietri 2002) shows that material with F^lOO would carry only a small 
fraction of the total energy, insufficient to account for the high observed 
values of E^^isQ. This is because the energy in the explosion is deposited into 
the bulk of the ejecta, and only a decreasing fraction of the total energy is 
transfered to a smaller fraction of the ejected mass that is accelerated to 
subsequently higher velocities as the supernova shock propagates into the 
steep density gradient at the outer edge of the progenitor star. For a jet the 
ejected mass can be much smaller than the accreted mass so that (F) ^ 1 
is possible, and a good fraction of the total energy in the outflow may be 
deposited in material with F > 100 that can power the GRB. Moreover, the 
energy requirements for a jet are much less severe, since the same observed 
value of E^^iso implies a much smaller beaming corrected radiated energy E^y^ 
so that typically only a small fraction (^ 10"'^ — 10~^) of the total available 
energy is required to power the jet. 

Finally, a somewhat more direct line of evidence in favor of a narrowly 
collimated outflow comes from achromatic breaks seen in the afterglow light 
curves of many GRBs (Fruchter et al. 1999, Kulkarni et al. 1999, Harrison 
et al. 1999, Halpern et al. 2000, Price et al. 2001, Sagar et al. 2001)0 In fact, 
such a "jet break" in the afterglow light curve was predicted before it was 
observed (Rhoads 1997, 1999, Sari et al. 1999). The exact cause of the jet 
break is discussed in detail in Sect. 111.81 A brief explanation follows. Due 
to relativistic beaming (and quasi-radial velocities of the emitting material) 
most of the observed emission comes from a visible region of angle ~ 1/F 
around our line of sight. As the jet decelerates by sweeping- up the external 
medium, the angular size of the visible region increases, until the edge of the 
jet becomes visible when F drops below the inverse of its half-opening angle, 
and from this point onward the observer "misses" the flux from outside the 

t Somewhat earlier, GRB 980519 showed a steep afterglow flux decay rate (~ t~^), which was 
interpreted as being due an earlier (unseen) jet break (Halpern et al. 1999) 
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edge of the jet, which would have been present if the flow were spherical, 
resulting in a faster flux decay (i.e., a jet break). 

For a bipolar jet that is uniform within its initial half-opening angle, ^O; 
the collimation or beaming factor is given by /b = 1 — cos^o ^ ^o/^' 
values of fh ~ 10^^ — 10^^ (corresponding to 6*0 values ranging between a 
few degrees and a few tens of degrees) have been derived in the pve-Swift 
era from such jet breaks (e.g., Frail et al. 2001). However, over the first 
years of Swift operations, only a handful of convincing jet-break candidates 
were identified (Blustin et al. 2006, Stanek et al. 2007, Willingale et al. 
2007, Kocevski &; Butler 2008), leading to concerns about the viability of 
this picture. Since then, however, deep optical and X-ray observations have 
revealed evidence for jet breaks in several additional Swift afterglows (Dai 
et al. 2008; Liang et al. 2008; Racusin et al. 2009; Burrows et al. 2009; Tanvir 
et al. 2010). A significant fraction of these aftreglows were monitored to late 
times, suggesting that jet breaks for typical Swift GRBs may be occurring 
at late times and faint flux levels that are beyond the limit of the standard 
ground- and space-based campaigns. This, in turn, might be attributed to 
Swiffs higher sensitivity compared to previous instruments, which results 
in the detection of fainter GRBs that (on average) correspond to wider jets. 



11.4 The jet structure 

The lack of well resolved GRB images makes it very difficult for the jet struc- 
ture to be deduced. The initial Lorentz factor during the prompt gamma-ray 
emission is very high, Fq ^ 100, and therefore we observe emission mainly 
from very small angles, 9 ^ Fq ^ < 10~^ rad, relative to our line of sight. This 
is a result of relativistic beaming (i.e., aberration of light; see Sect. I11.2p . 
For this reason, the prompt gamma-ray emission probes only a very small 
region of solid angle ~ ttFq , or a fraction ~ ~ lO"'^ - lO'^-^ of the 

total solid angle. Thus, the prompt emission provides no information about 
the ejecta that are moving in other directions (i.e., at angles ^ Fq ^ from 
our direction): the outflow could be spherical, or concentrated in a conical 
jet of half-opening angle > Tq ^ (provided that our line of sight is inside 
the cone, at an angle of ^Fq^ from its edge). 

During the afterglow, however, the Lorentz factor F of the emitting ma- 
terial decreases with time, since the afterglow shock decelerates as it sweeps 
up the external medium. This causes the beaming (or light aberration) ef- 
fects to gradually become less severe, as F becomes more moderate, so that 
we observe afterglow emission from a wider range of angles (of ^ F~^ from 
our line of sight, and not just from material moving almost directly towards 
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US, at angles ^Fq^, as in the prompt emission). This enables us to probe 
the jet structure over increasingly larger angular scales. 

For a jet with axial symmetry, its structure can be described by the 
distribution of its total energy (excluding rest energy) content per solid 
angle, and initial Lorentz factor Fq, as a function of the jet's polar angle 
9. We note that Fo(^) affects mainly the prompt gamma-ray emission and 
early afterglow, since it is largely "forgotten" after the local deceleration time 
tdec{(^) ~ -Rdec(^)/2cFQ(0), while £{9) affects also the late time afterglow 
emission. Here i?dec(^) is the local deceleration radius where most of the 
energy in the outflow at an angle 9 is transferred to the shocked external 
medium; for a power-law external density, pext{r) = Ar''', it is given by[3 
-Rdec(^) w [(3 - k)£{9)/Ac^Tl{9)]^/^^-^l The structure of GRB jets is very 
important for deducing their event rate and total energy, as well as for 
requirements on the jet production mechanisms. 

The most popular model for the structure of GRB jets is the uniform jet 
(UJ, or "top hat") model (Rhoads 1997, 1999, Panaitescu & Meszaros 1999, 
Sari et al. 1999, Kumar &: Panaitescu 2000, Moderski et al. 2000, Granot 
et al. 2001, 2002, Ramirez-Ruiz & Madau 2004, Ramirez-Ruiz et al. 2005), 
where £ and Fq are uniform within some finite half-opening angle, 9j, and 
sharply drop outside of 9y An alternative jet structure that is also rather 
popular is the universal structured jet (USJ) model (Meszaros et al. 1998, 
Lipunov et al. 2001, Rossi et al. 2002, Zhang & Meszaros 2002), where £ and 
Fq vary smoothly with 9, as a power law outside of some narrow core angle, 
typically with equal energy per decade in 9, £ cc 9^^. In the UJ model the 
different values of the jet break time, tj, in the afterglow light curve arise 
mainly due to different 9j (and to a lesser extent due to different ambient 
densities). In the USJ model, all GRB jets are assumed to be identical, and 
the different values of arise mainly due to different viewing angles, ^obsj 
from the jet axis. In fact, the expression for tj is similar to that for a uniform 
jet with £ — )• £{9 = ^obs) and 9j — >• ^obs- 

An alternative jet structure that has been proposed in the literature is 
one with a Gaussian angular profile (Zhang & Meszaros 2002, Kumar &: 
Granot 2003). It may be thought of as a more realistic version of a uni- 
form jet, where the edges are smooth rather than sharp. A Gaussian jet, 
£{9) oc exp(— 0^/20^), is approximately intermediate between the UJ and 
USJ models. However, it is closer to the UJ model than to the USJ model 
with £ oc 9~'^ in the sense that the energy in the wings of a Gaussian jet 



t We refer here to the angular structure, and ignore the radial structure of the outflow, which 

becomes unimportant once most of the energy is transfered to the shocked external medium. 
X This expression assumes a "thin" shell, as discussed in Sect. 111.51 
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is much smaller than in its core, whereas for a USJ with 6 oc wings 
there is equal energy per decade (in 6) in the wings, and therefore the wings 
contain more energy than the core (by about an order of magnitude). 

Another jet structure that is gradually receiving more attention is a two- 
component jet (Pedersen et al. 1998, Frail et al. 2000, Berger et al. 2003, 
Huang et al. 2004, Peng et al. 2005, Racusin et al. 2008) with a narrow 
uniform jet of initial Lorentz factor Pq ^ 100 surrounded by a wider uniform 
jet with Po ~ 10 — 30. Such a jet structure was predicted in the context 
of the cocoon in the collapsar model (Ramirez-Ruiz et al. 2002) and in 
the context of a hydromagnetically driven neutron-rich jet (Vlahakis et al. 
2003). This model has been invoked in order to account for sharp bumps 
(i.e., fast rebrightening episodes) in the afterglow light curves of GRB 030329 
(Berger et al. 2003) and XRF 030723 (Huang et al. 2004), but detailed cal- 
culations show that it cannot produce very sharp features in the light curve 
(Granot 2005). A different motivation for such jet structure comes about 
from the energetics of GRBs and X-ray flashes, which could help reduce the 
high efficiency requirements from the prompt gamma-ray emission (Peng 
et al. 2005). Later ^wi/t observations (e.g., Nousek et al. 2006) showed that 
while it can reproduce the early X-ray afterglow light curves, and specifically 
the shallow decay phase, it does not significantly help reduce the required 
gamma-ray efficiency (Granot et al. 2006). A different jet structure that has 
been suggested is that of a jet with a cross section in the shape of a "ring," 
sometimes referred to as a "hollow cone" (Levinson &; Eichler 1993, 2000, 
Eichler & Levinson 2003, 2004, Lazzati & Begelman 2005), which is uniform 
within Oc < < 6c + A6 where AO <^ 9c. 

Finally, it has been argued that there might be random variations on small 
angular scales in the prompt GRB brightness and in the energy per solid 
angle in the jet (and therefore also in the corresponding afterglow bright- 
ness), around some uniform mean value - the "patchy shell" model (Kumar 
&: Piran 2000). This model predicts wide variations in -E^^iso between dif- 
ferent lines of sight relative to the same GRB jet, as well as fluctuations in 
the afterglow light curve whose amplitude decreases with time as the visible 
region (of angle < P"^ around the line of sight) increases, thus effectively 
averaging the emission over an increasing number of bright and dim regions. 
The afterglow flux for different lines of sight approaches the same value at 
late times, thus implying a relatively narrow distribution of the afterglow 
(isotropic equivalent) luminosity compared to Ej^iso, and no obvious cor- 
relation with -E-Y,iso- While this seemed plausible at the time, later Swift 
observations showed a clear correlation where GRBs with large -E^^iso tend 
to have a more luminous afterglow emission (Nousek et al. 2006, Gehrels 
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et al. 2008). A more extreme version of the "patchy shell" model is the 
"mini-jets" model (Yamazaki et al. 2004), where the regions of bright emis- 
sion and high energy per solid angle are considered as discrete mini-jets, 
and the regions between them are assumed to have a negligible emission or 
energy per solid-angle. The main difference is that in this model many lines 
of sight fall between the mini-jets and the early emission is dominated by 
the one or more mini-jets closest to the line of sight, whose beaming cone 
initially does not encompass the line of sight. The original motivation for 
this model was unifying short GRBs with long GRBs and XRFs, where long 
GRBs, short GRBs and XRFs correspond to lines of sight with many, one 
and no mini-jets, respectively. However, the later discovery that the host 
galaxies of short GRBs are a different population from those of long GRBs 
(Gehrels et al. 2005, see also Fong et al. 2010) rules out the inclusion of 
short GRBs in such a unification scheme. 

Constraining the jet structure: efforts have been made to constrain the 
jet structure through statistical studies, linear polarization, afterglow light 
curves, off-axis viewing and orphan afterglows. Statistical studies have fo- 
cused on the log N — log S distribution where N is the number of GRBs 
above some threshold peak photon flux S (e.g., Firmani et al. 2004, Guetta 
et al. 2005), as weU as on dN/de (Perna et al. 2003) or dN/dzdO (Nakar 
et al. 2004), where z is the redshift and 9 is 6'j for the UJ model or ^obs 
for the USJ model, and is determined from the jet break time in the after- 
glow light curve. These studies were inconclusive, though they showed some 
preference for the UJ model over the USJ model. The same holds for stud- 
ies involving the shape of afterglow light curves (Granot & Kumar 2003), 
but nevertheless some jets structures could be ruled out using this method 
(Granot 2005). Studies based on afterglow polarization evolution need to 
deconvolve the effects of the jet structure from those of the magnetic field 
configuration in the post-shock emitting region, which make it difficult to 
draw strong conclusions about the jet structure without making similarly 
strong assumptions about the structure of the magnetic field itself. 

Orphan afterglows are events in which the late time afterglow is detected 
while the prompt gamma-ray emission is not. Unfortunately, no such events 
have been clearly observed so far, and the upper limits on the rates of 
orphan afterglows are still not very constraining for most of the different jet 
structures. Nevertheless, upper limits on the fraction of SN Ib/c that show 
late time radio emission of possible afterglow origin (Soderberg et al. 2006) 
already argue quite strongly against the presence of energetic, narrow jets 
(of half-opening angle 9q ~ I/Fq ~ 10~^'^) which would imply that a large 



12 



J. Granot & E. Ramirez-Ruiz 



fraction SNe Ib/c should harbor GRB jets (most of which point away from 

us and can be detected only at late times in the radio, when the jets become 
sub-relativistic). Such extremely narrow jets would also imply a very low 
true energy, which would be inconsistent with the lower limits on the kinetic 
energy inferred from late time radio afterglow observations (Berger, Kulkarni 
& Frail 2004, Frail et al. 2005) . On the other hand, two Type Ic supernovae 
have been recently observed to harbor relativistic jets (SN 2007gr - Paragi 
et al. 2010, and SN 2009bb - Soderberg et al. 2010), albeit with an energy 
either much (SN 2007gr) or slightly (SN 2009bb) lower than is needed to 
produce a typical hTight GRB. These two supernovae amount to about ~ 1% 
of relatively nearby Type Ib/c supernovae, for which a search for such jets 
has been performed at radio wavelengths. 

11.5 Dynamics of gamma-ray burst outflows 

Here we provide a brief summary of the key aspects for the dynamics of 
the interaction between the GRB outflow and the surrounding medium, and 
provide a simple intuitive explanation for some of the most important results. 
We begin with the relatively simpler case of a spherical outflow, and consider 
a uniform unmagnetized shell of ejecta of initial width Aq and Lorentz factor 
Tq in the lab frame, propagating into an external density pi{R)- A forward 
shock is driven into the ambient medium, while the ejecta are decelerated by 
a reverse shock. Thus four regions exist: (1) unshocked external medium, 
(2) shocked external medium, (3) shocked ejecta, and (4) unshocked (freely 
expanding) ejecta. All the velocities are measured relative to region 1, while 
the pressure p and rest mass density p (or number density n) are measured 
in the fluid rest frame. A subscript i between 1 and 4 refers to region i, 
while a subscript ij refers to the relative velocity of regions i and j, so that 
Til = Tii = Fj. Given F4 = Fq, Fi = 1, p4, pi, and assuming the shell of 
ejecta and external medium are both cold {pi <^ pi(? for i = 1, 4), there are 
8 unknown hydrodynamic quantities (p2, P2, ^2, Ps, P3, Ta, and the Lorentz 
factors of the forward and reverse shock fronts) that can be found from 
the conditions across the contact discontinuity separating regions 2 and 3 
{P2 = P3 and F2 = F3 = F) as well as the shock jump conditions (continuity 
of the energy, momentum, and particle fluxes) across the forward and reverse 
shocks (between regions 1 & 2, and 4 &; 3, respectively). For simplicity, this 
is treated in planar symmetry, and the spherical nature of the flow enters 
only when the evolution of the flow with radius is considered. 

An approximate, order of magnitude, estimate of the dynamics at this 
stage can be obtained by equating the ram pressure of the incoming fluid 
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from regions 4 and 1, as seen from the contact discontinuity (the rest frame 
of regions 2 and 3): piuli ~ P4U34, where u = VP is the 4- velocity. For 
P4 = pi we must have U21 = U34 due to symmetry, and for r4 = To 3> 1 
this imphes that both the forward and reverse shocks are relativistic, U21 = 
U43 = [{^4 - l)/2]^/^ « (r4/2)^/2 » 1. As long as the forward shock is 
relativistic, 1 <C = r2 ~ U21 ~ U34{p4/ pi)^^^; when the reverse shock is 
relativistic then 1 < U43 r4/2r3, so that r2 « (ro/2)^/2(p^/p^)i/4 

and ~ (ro/2)^/^(pi/p4)^/'^ (Sari &: Piran 1995). Therefore, the condition 
for the forward shock to be relativistic is Fq ^ Pi/P4, which is typically 
always satisfied (until very late times when the flow becomes Newtonian), so 
that U21 ~ F2 and 1*43 ~ T2{pi/ Pa)^^"^ ■ The condition for the reverse shock to 
be relativistic is Fq ^ pi/ pi- The revere shock is Newtonian in the opposite 
limit, Fq <C P4/P1, and, in this case, F2 ~ Fq so that M43 ~ Fo(pi/y04)^/^ <C 1. 

Note that U43 m.m{a^/^,a^/^) where a = F§pi/p4 - Ac'^T^R'^-''A/E ~ 
FqE!^^^ A/l^~^ for a power-law external density, pi = AR~^, where I ~ 
(£;/^c2)V(3-fc) is the Sedov length. For a narrow distribution of Lorentz 
factors A pa Aq = const, and therefore a oc i?^"*^, i.e., for A; < 2 the 
reverse shock is initially Newtonian and strengthens with radius. If there is 
a reasonable spread in the Lorentz factor of the outflow, AFq ~ Fq, then the 
shell can spread, A ~ max(Ao, R/Tq), where A (x R and therefore a oc 
at R > Rg ^ TqAo so that then the reverse shock strengthens for A; < 3. 

If the reverse shock is relativistic by the time it finishes crossing the shell, 
then most of the energy is transfered to the shocked external medium within 
a single shell crossing. If there is only a very small spread in Fq then the 
reverse shock can still be Newtonian when it finishes crossing the shell. In 
this case a large number of Newtonian shocks and rarefaction waves may 
need to cross the shell before most of the energy is transfered to the shocked 
external medium (Sari 1997). However, if there is a reasonable spread in 
the Lorentz factor of the outflow, AFq ~ Fq, then the shell starts to spread 
before the reverse shock finishes crossing it, in such a way that by the time it 
crosses the shell the reverse shock already becomes mildly relativistic, so that 
most of the energy is transfered to the shocked external medium in a single 
shell crossing time (Sari &; Piran 1995). The dividing line between these 
two cases corresponds to 1 = a{Rs) ~ Fq*-"^ ''\Aq/1)^~'^, where a{Rs) > 1 
implies a relativistic reverse shock or a "thick" shell and a{Rs) < 1 implies 
a Newtonian reverse shock (without spreading) or a "thin" shell. 

Most of the energy is transferred to the shocked external medium at a 
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radius i?dec ~ J'^dec^^~'''' where T^ec = ^{Rdec) ~ min[ro,rcr] and 
, x(3-/c)/2(4-fc) (28Qe/'Eli''no'^%-J'^' {k = 0), 

rer=(^) =<! 



A, 



(11.3) 

i 70CV^£;V4a;VV-;/^ (fc = 2) , 

where C = (1 + ^)/2, Tgrb = (1 + z)Ao/c = SOTsoS is the duration of the 
GRB, E = lO^'^E'saerg is the (isotropic equivalent) energy of the flow, and 
n = no cm~^ = A/nip for A; = 0, while = A/ {5 x 10^"*^ gr cm~^) for k = 2. 
For To > r cr we have a thick shell or relativistic reverse shock, and the 
observed deceleration time is similar to the duration of the GRB, t^ec ~ 
Rdec/2cT'^^^ ~ (1 + z)Aq/c ^ Tgrb- For Fq < Fcr we have a Newtonian (or 
at most mildly relativistic) reverse shock or a thin shell, and in this case 

^dec ^ (Vc)ro^^^ '^^^^^ > Tgrb (where l/c ^ ^nr is the non-relativistic 
transition time) and is given by 

idec = (l + ^)^= (11.4) 
[ 5.%E53A-\To/m-U {k = 2). 

Once most of the energy in the GRB outflow is transfered to the shocked 
external medium, the flow becomes self-similar, and is described by the 
Blandford & McKee (1976) solution. The scaling of the basic quantities with 
radius during this stage may be understood as follows. In this stage most of 
the energy is in the shocked external medium, and therefore the dynamics 
may be found by energy conservation within this region. In the rest frame 
of the the shocked external medium (region 2), the cold upstream external 
medium approaches at a Lorentz factor = which for simplicity we 
denote here by F, and the velocities of the particles are randomized at the 
shock font, such that the ordered bulk motion in the upstream region is 
converted to random motion of the particles in the downstream region (2), 
with the same average Lorentz factor. Therefore, in the rest frame of the 
shocked downstream medium, the average energy per particle (including its 
rest energy) is F times its rest energy. In the lab frame the energy is larger 
by a factor of F, implying that the total kinetic energy in the lab frame is 
E = {T'^ - l)M{R)c^ = v?M{R)c^, where M{R) is the total swept-up rest 
mass up to radius R, u = TP is the 4- velocity, and we have deducted the rest 
energy in order to obtain the kinetic energy (only the latter is conserved, 
since new rest-mass is added to the shocked region as the external medium 
is swept-up by the forward shock). 

For an external density that varies as a power law with radius, pi oc R~^, 
we have E oc u^R^~^ and therefore conservation of energy in the lab frame 
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implies that u oc R^'^^^y^. For a relativistic flow u T and T oc i?('^~^)/2 
Eventually the flow becomes Newtonian and approaches the Sedov- Taylor 
solution. In this regime u ~ /3 and therefore /3 oc i?^'^'^'^)/^. Both scalings 
given above apply to the adiabatic case where energy losses (e.g., due to 
radiation) or gains (e.g., due to late time energy injection from the central 
source) can be neglected. For simplicity we consider here only < 3 for 
which u decreases with R. The pressure in the shocked region scales as 
p ~ pic^u^ oc R^"' , while the particle rest mass density (in the comoving 
frame) scales as p2 ~ Fpi, i.e., as in the relativistic case and as 

R~^ for the Newtonian case. A spherical flow becomes non-relativistic once 
r2 f« E/M{R)c^ ~ 1, at 

i/(3-fc) f 5.3C-E5^^no^^^yr (A; = 0) , 



I 



(3- fc)^ 



STC-EssA-iyr (^ = 2) 



(11.5) 



Variations on the relatively simple dynamics described above after the 
deceleration epoch involve relaxing either the assumption of a constant en- 
ergy or of spherical symmetry. The energy can decrease due to radiative 
losses or increase due to energy injection, either from a relativistic wind 
caused by prolonged activity of the central source, or slow material that was 
ejected promptly from the source but catches up with the shocked exter- 
nal medium as the latter decelerates to a somewhat smaller Lorentz factor 
(e.g., Blandford & McKee 1976, Cohen & Piran 1999, Sari & Meszaros 2000, 
Ramirez-Ruiz et al. 2001). Here we discuss deviation from spherical sym- 
metry, for which the dynamics are more complicated and less certain. For 
simplicity, we shall consider an axisymmetric double-sided jet that is uni- 
form within an initial half-opening angle of 6*0 <C 1 around its symmetry 
axis with sharp edges. 

If a signal (e.g., a sound or rarefaction wave) can travel laterally within 
the jet at a speed of /SgC in its comoving frame, then it traverses an angle 
of d9 « dR±/R = /3scdt'/R ^ MR/R^, so that even for ~ 1 (which 
may indeed be expected) information traverses an angle of ~ l/F on the 
dynamical (or radius doubling) time. Therefore, the center of the jet knows 
about its edge only when F drops to ~ 1/Oq; at earlier times it behaves 
as if it were part of a spherical flow with the same external density profile 
and isotropic equivalent energy Eiso- Using the spherical adiabatic scaling 
F oc i?('=-3)/2 a signal starting at Rq traverses A6' Ri 2/3s/(3 - k)r{R) by 
a radius R ^ Rq for k < 3. Therefore, when A9 = 6q corresponding to 
F^o ~ 2/3s/(3 — k), the center of the jet comes into causal contact with its 
edge, and the jet can in principle start to spread sideways rapidly. 
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Early analytic works (Rhoads 1997, 1999, Sari et al. 1999) assumed that 
the jet half-opening angle grows as = + ~ + /^s/r and indeed 
quickly expands laterally once V drops below ~ /3s/^o ~ l/^O; around the 
same time as the edges of the jet become visible^. As a result, a steepening in 
the afterglow light curve is produced, which is known as a jet break and was 
detected soon after it was predicted. The jet break occurs at a radius i?j ~ 

i?NR(-E') ~ ^0^^^ '^^i^NRl-E'iso), and the simple analytic models (e.g., Rhoads 
1999) suggest that at i? > T ~ O^^e^-^'^^ and 6'j ~ 9o{R^/ R)e^-^/^i 
until r ~ 1 at -R ~ i?j(l — In^o)- Different analytic models vary in their exact 
assumptions and results (e.g., Piran 2000, Granot 2007), but the exponential 
decrease in F and increase in 9^ with radius are a generic prediction. In this 
picture the jet becomes sub-relativistic at close to R^^{E) (to within a 
logarithmic factor), and therefore the transition to the late time asymptotic 
spherical Newtonian Sedov- Taylor solution is relatively prompt and smooth 
in this scenario. 

However, numerical studies (Granot et al. 2001, Kumar & Granot 2003, 
Cannizzo et al. 2004, Zhang &: MacFadyen 2009), and in particular full spe- 
cial relativistic two dimensional hydrodynamic simulations, which are the 
best and most reliable calculations to date, show that the lateral expansion 
of the jet is very modest as long as it remains relativistic. The big difference 
compared to the results of simple analytic models may be attributed to the 
over-simplified assumptions, such as the jet being perfectly uniform with 
sharp edges within some finite half-opening angle with a purely radial veloc- 
ity, which are shown to be invalid in the numerical simulations. Because of 
the lack of significant sideways expansion during the relativistic stage, the 
gross properties of the jet may be approximated to zeroth order by a conical 
section of half opening angle ~ 6q out of a spherical flow. Therefore, the jet 
becomes non-relativistic at a radius ~ i?NR(-E'iso) ~ 9^'^^^^ ''^Rtqji{E), i.e., 
near the Sedov length corresponding to its initial isotropic equivalent energy 
E'iso, which for a narrow jet (^o ^ 1) is significantly larger than the Sedov 
length corresponding to its true energy E Siso^o/^- This causes the tran- 
sition into a spherical flow described by the Sedov- Taylor solution (with the 
true energy E) to extend over a large range of times, with a rather modest 
growth in the (maximal) radius during this transition period (Granot et al. 
2005). 



f Since the visible region is within an angle of ~ l/F around the line of sight due to relativistic 
beaming and light travel effects. 
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11.6 The afterglow emission 

The dominant emission mechanism during the afterglow stage is believed to 
be synchrotron radiation, produced as relativistic electrons accelerated by 
the afterglow shock gyrate in the magnetic fields within the shocked external 
medium. A synchrotron origin of the afterglow emission is supported by 
the detection of linear polarization at the level of ~ 1% — 3% in several 
optical or NIR afterglows (see Sect. Ill.Op . and by the shape of the broad 
band spectrum, which consists of several power-law segments that smoothly 
join at some typical break frequencies (Galama et al. 1998). Synchrotron 
self-Compton (SSC) - the inverse-Compton scattering of the synchrotron 
photons to (much) higher energies by the same population of relativistic 
electrons that emits the synchrotron photons - can sometimes dominate the 
afterglow flux in the X-rays (Sari & Esin 2001, Harrison et al. 2001), and 
may affect the synchrotron emission by increasing the electron cooling. 

Since the physics of relativistic collisionless shocks (and in particular how 
they amplify the magnetic field and accelerate particles to a non-thermal 
relativistic energy distribution) are still not well understood from first prin- 
ciples, simple assumptions are usually made that conveniently parameterize 
our ignorance. The electrons are usually assumed to be (practically in- 
stantaneously) shock-accelerated into a power-law distribution of energies, 
dN/dje oc for > Im, and thereafter to cool both adiabatically and 
due to radiative losses. It is further assumed that practically all of the elec- 
trons take part in this acceleration process and form such a non-thermal 
(power-law) distribution, leaving no thermal component (which is not at all 
clear or justified; e.g., Eichler & Waxman 2005). The relativistic electrons 
are assumed to hold a fraction eg of the internal energy immediately behind 
the shock, while the magnetic field is assumed to hold a fraction of the 
internal energy everywhere in the shocked region. 

The spectral emissivity in the comoving frame of the emitting shocked 
material is typically approximated as a broken power law (in some cases 
the more accurate functional form of the synchrotron emission is used, e.g., 
Wijers & Galama 1999, Granot &: Sari 2002). Most calculations of the 
light curve assume emission from an infinitely thin shell, which represents 
the shock front (some integrate over the volume of the shocked fiuid taking 
into account the appropriate radial profile of the flow, e.g., Granot &: Sari 
2002; see Figure [Il.3|) . One also needs to account for the different arrival 
times of photons to the observer from emission at different lab frame times 
and locations relative to the line of sight, as well as the relevant Lorentz 
transformations of the emission into the observer frame. SSC is included 
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Fig. 11.3. The afterglow synchrotron spectrum, calculated for the Blandford & Mc- 
Kee (1976) spherical self-similar solution, under standard assumptions, using the 
accurate form of the synchrotron spectral emissivity and integration over the emis- 
sion from the whole volume of shocked material behind the forward (afterglow) 
shock (for details see Granot & Sari 2002). The different panels show the five pos- 
sible broad band spectra of the afterglow synchrotron emission, each corresponding 
to a different ordering of the spectral break frequencies. Each spectrum consists of 
several power-law segments (PLSs; each shown with a different color and labeled 
by a different letter A-H) that smoothly join at the break frequencies (numbered 
1-11). The broken power-law spectrum, which consists of the asymptotic PLSs 
that abruptly join at the break frequencies (and is widely used in the literature), 
is shown for comparison. Most PLSs appear in more than one of the five different 
broad band spectra. Indicated next to the arrows are the temporal scaling of the 
break frequencies and the flux density at the different PLSs, for a uniform (ISM; 
fc = 0) and stellar wind (WIND; k — 2) external density profile. 
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in some (but not all) works, even though it can also effect the synchrotron 
emission through the enhanced radiative cooling of the electrons. 



11.7 Afterglow images 

In order to calculate the afterglow image, or how the afterglow would appear 
on the plane of the sky if it were angularly resolved, we need to specify the 
dynamics, in addition to an emission model, here assumed to be synchrotron 
radiation. As discussed in Sect. 111.51 at early times before the jet break time 
in the afterglow light curve the dynamics may be reasonably approximated 
as part of a spherical flow that is described by the self-similar Blandford & 
McKee (1976) solution. Here we shall concentrate on this stage, for which 
the afterglow image is also self-similar. 

During the self-similar spherical evolution stage (before the jet break time, 
for a jet), the afterglow image has circular symmetry around the line of sight 
(where the surface brightness depends only on the distance from the center 
of the image) , and is confined within a circle on the sky with a radius 



3-91 (ff) (fe) (fc = 0) 



IQi'^cm I /p \l/4/t^ \3/4 

2.39 (f^) (fe) {k 



(11.6) 



(see Figure [TO|) . 

This corresponds to an angular radius of 



(11.7) 



where dA{z) = 10^''-'^(iA,27.7 cm^ 

More generally, the afterglow image size during the self-similar spheri- 
cal stage scales with the observed time as oc The image 
size grows super-luminally with an apparent expansion velocity of T^-^{R^)c. 
The expected afterglow images during this self-similar regime are shown in 
Figures 111.51 and 111.61 The normalized surface brightness profile within the 
afterglow image is independent of time due to the self-similar dynamics, and 
changes only between the different power-law segments of the synchrotron 
spectrum, and for different external density profiles. The image becomes 
increasingly limb-brightened at higher frequencies, and for smaller values of 
k (Granot & Loeb 2001, Granot 2008). 

f For a standard cosmology (fijvf = 0.27, Qa = 0.73, h = 0.71) dj\{z) has a maximum value of 
5.45 X 10^'' cm {dA,27.7 = 1-09) for z = 1.64. 
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Fig. 11.4. Schematic illustration of the equal arrival time surface (thick black line), 
namely, the surface from where the photons emitted at the shock front arrive at 
the same time at the observer (on the far right-hand side). The maximal lateral 
extent of the observed image, i?^, is located at an angle , where the shock radius 
and Lorentz factor are i?* and = Fsh(-R*), respectively. The area of the image 
on the plane of the sky is S±^ = itR\. The shock Lorentz factor, Fgh, varies with 
radius R and angle 6 from the line of sight along the equal arrival time surface. 
The maximal radius Ri on the equal arrival time surface is located along the line 
of sight. If, as expected, Fgh decreases with i?, then F/ = Tsh{Ri) is the minimal 
shock Lorentz factor on the equal arrival time surface, (from Granot et al. 2005). 

Below the self-absorption frequency, the specific intensity (surface bright- 
ness) represents the Rayleigh-Jeans portion of a black-body spectrum with 
the blue-shifted effective temperature of the electrons at the corresponding 
radius along the front side of the equal arrival time surface of photons to the 
observer {R^ < R < Ri m Figure III .41) . Above the cooling break frequency 
the emission originates from a very thin layer behind the shock front, where 
the electrons whose typical synchrotron frequency is close to the observed 
frequency have not yet had enough time to significantly cool due to radia- 
tive losses. This results in a divergence of the surface brightness at the outer 
edge of the image (Sari 1998, Granot & Loeb 2001). 

After the jet break time the afterglow image is no longer symmetric around 
the line of sight to the central source for a general viewing angle (which is 
not exactly along the jet symmetry axis), and its details depend on the 
hydrodynamic evolution of the jet (so that in principle it could be used 
to constrain the jet dynamics). Therefore, a realistic calculation of the 
afterglow image during the more complicated post-jet break stage requires 
the use of hydrodynamic simulations, and still remains to be done. 

The afterglow image may be indirectly resolved via gravitational lensing 




Fig. 11.5. The afterglow images for different power-law segments of the spectrum, 
for a uniform (fc = 0) and wind (k = 2) external density profile (from Granot & Loeb 
2001), calculated for the Blandford & McKee (1976) spherical self similar solution, 
using the formalism of Granot & Sari (2002). Shown is the surface brightness, 
normalized by its average value, as a function of the normalized distance from the 
center of the image, r ~ RsinO/ R± (where r = at the center and r = 1 at the 
outer edge). The image profile changes considerably between different power-law 
segments of the afterglow spectrum, Fi, oc . There is also a strong dependence 
on the density profile of the external medium, pext oc R~^. 
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Fig. 11.6. An illustration of the expected afterglow image on the plane of the sky, 
for three different power-law segments of the spectrum (from Granot et al. 1999a, b), 
assuming a uniform external density and the Blandford McKee (1976) self-similar 
solution. The image is more limb brightened at power-law segments that correspond 
to higher frequencies. 

by a star in an intervening galaxy (along, or close to, our line of sight to 
the source). This is because the angular size of the Einstein radius (i.e., the 
region of large magnification around the lensing star) for a typical star at a 
cosmological distance is ~ 1 /ias (hence the name micro-lensing) - compara- 
ble to the afterglow image size after a day or so. Since the afterglow image 
size grows very rapidly with time, different parts of the image sample the 
regions of large magnification (close to the point of infinite magnification 
just behind the lensing star) with time, and therefore the overall magnifica- 
tion of the afterglow flux as a function of time probes the surface brightness 
profile of the afterglow image. This results in a bump in the afterglow light 
curve, which peaks when the limb-brightened outer part of the image sweeps 
past the lensing star; the peak of the bump is sharper when the afterglow 
image is more limb-brightened (Granot & Loeb 2001). It has been suggested 
that an achromatic bump in the afterglow light curve of GRB 000301C af- 
ter ~ 4 days might have been due to micro-lensing (Garnavich et al. 2000). 
If this interpretation is true, then the shape of the bump in the afterglow 
light curve requires a limb-brightened afterglow image, in agreement with 
theoretical expectations (Gaudi et al. 2001). 



11.8 The cause of the jet break 

The jet break in the afterglow light curve has been argued to be the com- 
bination of (i) the edge of the jet becoming visible, and (ii) fast lateral 
spreading. Both effects are expected to take place around the same time. 
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when the Lorentz factor, F, of the jet drops below the inverse of its initial 
half-opening angle, ^o- This can be understood as follows. 

When r drops below Oq^, the edge of the jet becomes visible. This is be- 
cause relativistic beaming limits the region from which a significant fraction 
of the emitted radiation reaches the observer to within an angle of ~ 
around the line of sight (0^r^^). Once the edge of the jet becomes visi- 
ble, then, if there is no significant lateral spreading, only a small fraction, 
{Ve^f < 1, of the visible region is occupied by the jet and, as a result, 
there would be "missing" contributions to the observed flux compared to 
a spherical flow. This would cause a steepening in the light curve, i.e., 
a jet break, where the temporal decay index asymptotically increases by 
Aa = (3 — k)/{A — k) (since the fraction of the visible region occupied by 
the jet is (r6'j)2 ~ 

When r drops below Oq^ , the center of the jet comes into causal contact 
with its edge, and the jet can in principle start to expand sideways signifi- 
cantly. It has been argued that at this stage it would indeed start to expand 
sideways rapidly, at close to the speed of light in its own rest frame. In this 
case, during the rapid lateral expansion phase the jet opening angle grows 
as ~ and exponentially with radius (see Sect. [TT3]) . This causes the 
energy per unit solid angle, in the jet to drop with observed time, and the 
Lorentz factor to decrease faster as a function of observed time; this gives 
rise to a steepening in the afterglow light curve compared to a spherical flow 
(for which £ remains constant and F decreases more slowly with the observed 
time). However, in this large fraction of the visible region remains 

occupied by the jet (since TO-^ remains ~ 1), so that the first cause for the jet 
break (the edge of the jet becoming visible, and the "missing" contributions 
from outside the edge of the jet) is no longer important. Therefore, for fast 
lateral spreading, the jet break occurs predominantly both as a result of 
the energy per solid angle £ decreasing with time and the Lorentz factor 
decreasing with observed time faster than for a spherical flow. 

It is important to keep in mind, however, that numerical studies show that 
the lateral spreading of the jet is very modest as long as it is relativistic (see 
Sect. 111.5]) . This implies that lateral spreading cannot play an important role 
in the jet break, and the predominant cause of the jet break is the "missing" 
contribution from outside of the jet, once its edge becomes visible. 

A potential problem with this description is that if the jet half-opening 
angle remains roughly constant, 6^ ~ ^o? the asymptotic change in the tem- 
poral decay index is only Aa = 3/4 for a uniform external medium {k = 0) 
or even smaller for a wind (Aa = 1/2 for k = 2), while the values inferred 
from observations are in most cases larger (see Figure 3 of Zeh et al. 2006). 
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Fig. 11.7. The temporal decay index a as a function of the observed time (in days) 
across the jet break in the hght curve, for p = 2.5. Top panel: results in the spec- 
tral range v > max(:^m,'^c) using a semi-analytic model with no lateral spreading 
(Granot 2005), for a uniform (fc = 0, next — Icm^"^) and wind (fc = 2, A, = 1) 
external density profile, with 0q = 0.1 and i?k,iso — 2 x 10^^ erg. Bottom panel: re- 
sults for the spectral range i^m < v < i^c, for = 0-2 and a uniform density (fc — 0, 
riext = 1 cm~^, £^k,iso = 10^^ erg); the figure compares the result of a semi-analytic 
model (Granot 2005) to those of a hydrodynamic simulation (Granot et al. 2001). 
In both panels the dashed lines show the asymptotic values of a before and after the 
jet break, for a uniform jet with no lateral spreading, for which Aa = (3 — fc)/(4— fc). 



This apparent discrepancy may be reconciled as follows. While the asymp- 
totic steepening is indeed Aa = (3 — A;)/(4 — k) when lateral expansion is 
negligible, the value of the temporal decay index a (where F^^ oc t~") ini- 
tially overshoots its asymptotic value. Since the temporal baseline that is 
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used to estimate the post-jet break temporal decay index 02 is typically no 
more than a factor of several in time after the jet break tim^, ij, the value 
of a during this time is larger than its asymptotic value, 02- This causes 
the value of Aa that is inferred from observations to be larger than its true 
asymptotic value. 

The overshoot in the value of a just after the jet break time can be seen 
in Figure 111.71 and is much more pronounced in the light curves calcu- 
lated using the jet dynamics from a hydrodynamic simulation, compared to 
the results of a simple semi-analytic model. This overshoot occurs because 
the afterglow image is limb-brightened (see Figure Ill.Sp and therefore its 
brightest part - its outer edges - is the first region whose contribution to the 
observed flux is "missed" as the edge of the jet becomes visible. The over- 
shoot increases for increasingly limb-brightened afterglow images (e.g.,, for 
u > max[i^m, Vc\ in the upper panel of Figure [TL71 compared io Vm < v < Vc 
in the lower panel of Figure [11.7p . For a wind density {k = 2) the limb- 
brightening is smaller compared to a uniform density {k = 0), at the same 
power-law segment of the spectrum (see Figure [Il.5|) . and the Lorentz fac- 
tor r decreases more slowly with the observed time. Because of this no 
overshoot is seen in the semi-analytic model shown in the upper panel of 
Figure [11. 71 for a wind density profile (fc = 2), and the jet break is smoother 
and extends over a larger factor in time. The asymptotic post-jet break 
value of the temporal decay index (02) is approached only when the part 
of the image of the corresponding spherical flow that is occupied by the 
jet covers only the relatively uniform central part of the image, and not its 
brighter outer edge. 

The jet break in light curves calculated from hydrodynamic simulations is 
sharper than in semi-analytic models (where the emission is taken to be from 
a 2D surface - usually a section of a sphere within a cone). In semi-analytic 
models the jet break is sharpest with no lateral expansion, and becomes more 
gradual with faster assumed lateral expansion. For example, in the lower 
panel of Figure lll.7[ where the viewing angle is along the jet axis and the 
external density is uniform, most of the change in the temporal decay index 
a occurs over a factor of 2 in time for the numerical simulation, and over a 
factor of ~ 3 in time for the semi-analytic model (which assumes no lateral 
expansion; the jet break would be more gradual with lateral expansion). 
For both types of models, the jet break is more gradual and occurs at a 
somewhat later time for viewing angles further away from the jet symmetry 
axis but still within its initial opening angle, although this effect is somewhat 



t This is usually because the flux becomes too dim to detect above the host galaxy, or since a 
supernova component becomes dominant in the optical. 
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more pronounced in semi-analytic models (Granot et al. 2001, Rossi et al. 
2004). 

11.9 Polarization: afterglow and prompt 

Afterglow: Linear polarization at the level of a few percent has been de- 
tected in the optical and NIR afterglow of several GRBs (see Covino et al. 
2004 and references therein). This was considered to be a confirmation that 
synchrotron radiation is the dominant emission mechanism in the afterglow. 
The synchrotron emission from a fluid element with a locally uniform mag- 
netic field is linearly polarized in the direction perpendicular to the projec- 
tion of the magnetic field onto the plane normal to the wave vector. Since 
the source moves relativistically, aberration of light (see Sect. 111.^ must be 
accounted for when calculating the observed local direction of polarization. 
Figure [T 1 . 81 shows the predicted local polarization map from emission by an 
ultra-relativistic expanding shell, for two different magnetic field structures: 
a magnetic field that is random within the plane normal to the radial di- 
rection {left panel) as could be expected from a shock produced magnetic 
field (e.g., Medvedev & Loeb 1999), and an ordered magnetic field normal 
to the radial direction (as could be expected in the prompt emission for a 
magnetic field coherent on angular scales ^> I/Tq that is advected from the 
central source). 

Since the afterglow image is almost always never resolved, we can only 
measure the average polarization over the whole image. For this reason, a 
shock produced magnetic field that is symmetric about the shock normal 
will procure no net polarization for a spherical fiow (since the polarization 
pattern across the image is in this case symmetric around the center of the 
image, and the polarization averages out to zero when summing over the 
contributions from the whole image). For a shock-produced magnetic field 
one thus needs to break the spherical symmetry of the flow to produce net 
polarization. The simplest and most natural way of doing this is considering 
a jet, or narrowly collimated outflow (e.g.. Sari 1999, Ghisellini & Lazzati 
1999). In this picture a jet geometry together with a line of sight that is 
not along the jet symmetry axis (but still within the jet aperture, in order 
to see the prompt GRB) is needed to break the symmetry of the afterglow 
image around our line of sight. 

For a uniform jet (the UJ model) this predicts two peaks in the polar- 
ization light curve around the jet break time tj, if F^j < 1 decreases with 
time at t > ij (Ghisellini & Lazzati 1999, Rossi et al. 2004), or even three 
peaks if F^j ~ 1 at i > t j (Sari 1999), where in both cases the polariza- 
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Fig. 11.8. The predicted polarization map for synchrotron emission from a thin 
spherical ultra-relativistic shell expanding with a Lorentz factor F ^ 1 . The double- 
sided arrows show the direction of the linear polarization (the wave electric vector), 
while their length depends monotonically on the polarized intensity (in a non-trivial 
way, for display purposes). The circle indicates an angle of 1/F around the line of 
sight, which is the region responsible for most of the observed flux. The left panel 
is for a magnetic field that is random within the plane of the shell (normal to the 
radial direction), for which the polarization direction always points at the center of 
the image (corresponding to the line of sight to the center of the spherical shell), 
where the polarization vanishes (due to symmetry consideration) . The right panel 
is for an ordered magnetic field within the plane of the shell that is coherent over 
angular scales » 1/F (Granot & Konigl 2003). In this case the direction of the 
ordered magnetic field clearly breaks the symmetry around the center of the image, 
resulting in a large net polarization. For simplicity, the map is for a constant 
emission radius, rather than for a constant photon arrival time. 

tion vanishes and reappears rotated by 90° between adjacent peaks. The 
latter is a distinct signature of this model. For a structured jet (the US J 
model), the polarization position angle is expected to remain constant in 
time, vi^hile the degree of polarization peaks near the jet break time tj (Rossi 
et al. 2004). A similar qualitative behavior is also expected for a Gaussian 
jet, or other jet structures with a bright core and dimmer wings (although 
there are obviously some quantitative differences). 

The different predictions for the afterglow polarization light curves for dif- 
ferent jet structures imply that afterglow polarization observations could in 
principle help constrain the jet structure. In practice, however, the situation 
is much more complicated, mainly since the observed polarization depends 
not only on the jet geometry, but also on the magnetic field configuration 
within the emitting region, which is not well known. For example, an or- 
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dered magnetic field component in the emitting region (e.g., due to a small 
ordered magnetic field in the external medium) may dominate the polarized 
flux, and therefore the polarization light curves, even if it is sub-dominant 
in the emitting region compared to a random (shock generated) magnetic 
field component in terms of the total energy in the magnetic field (Granot 
&; Konigl 2003). Other models for afterglow polarization include a magnetic 
field that is coherent over patches of a size comparable to that of causally 
connected regions (Gruzinov & Waxman 1999), and polarization that is in- 
duced by microlensing (Loeb & Perna 1998) or by scintillations in the radio 
(Medvedev & Loeb 1999). 

The last two models involve mechanisms that give a different weight to 
emission from different parts of the afterglow image, and thus break its 
symmetry around the line of sight for a shock produced magnetic field (so 
that the polarization from the different parts of the image no longer fully 
cancels out). While microlensing or radio scintillations are external to the 
source, there are also mechanisms intrinsic to the source that can produce 
a similar effect on the polarization. An example is a "patchy shell" , where 
the energy per solid angle and the corresponding afterglow brightness vary 
randomly with the lateral location within the flow (see Sect. I11.4p . As 
the afterglow shock decelerates and the angular size of the visible region (of 
< around the line of sight) increases, new bright and dim regions become 
visible and cause correlated variability in the total afterglow flux (bumps or 
wiggles in the afterglow light curve) and in its linear polarization - both 
in the degree of polarization and its position angle (Nakar & Oren 2004). 
Soon after this was predicted (Granot & Konigl 2003)0 such a correlated 
variability in the afterglow light curve and its polarization was reported for 
two particularly variable optical afterglows, of GRB 021004 (Rol et al. 2003) 
and GRB 030329 (Greiner et al. 2003). 

Prompt Emission: in the prompt soft gamma-ray emission the nature 
of the dominant emission mechanism remains uncertain. Moreover, it is very 
difficult to measure the polarization at such photon energies (hard X-rays or 
soft gamma-rays), compared to the optical, NIR, or radio. There have been 
some claims (Coburn & Boggs 2003, Wilhs et al. 2005, McGlynn et al. 2007, 
Gotz et al. 2009) of detection of a high degree of linear polarization in the 
prompt gamma-ray emission of some GRBs (with a particularly high fluence, 
as good photon statistics are crucial for such measurements). However, these 
claims remain rather controversial (Rutledge & Fox 2004, Wigger et al. 



t A correlated variability in the afterglow light curve and its polarization is also expected in other 
scenarios, such as microlensing, radio scintillation and non-axisymmetric "refreshed shocks" or 
variations in the external density. 
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2004). A detection of polarization during the prompt emission phase could 
nonetheless help constrain the dominant emission mechanism. 

An ordered magnetic field within the outflow on angular scales ^ l/F 
can produce a high net polarization, of tens of percent (Granot & Konigl 
2003, Granot 2003, Lyutikov et al. 2003), which is only slightly smaller 
than the maximal polarization of the local synchrotron emission, as the 
polarization vector is more or less aligned in the regions that contribute 
the most to the observed flux (see right panel of Figure Ill.Sp . Such a high 
polarization can also be produced by synchrotron emission from a shock 
generated magnetic field (Waxman 2003, Granot 2003, Nakar et al. 2003) 
or bulk inverse-Compton scattering of an external photon field (Shaviv & 
Dar 1995, Lazzati et al. 2004). However, these mechanisms require the line 
of sight to be near the edge of the jet (to within an angle of < 1/r). Q In 
order for such a line of sight not to be too rare, a narrow jet is required: 
^o^a few r^^. This requirement may be viable for the brightest GRBs, 
for which the prompt polarization can be measured, as such bright events 
usually correspond to very narrow jets. Nevertheless, statistical studies over 
a sample of GRBs, or time resolved polarimetry of different emission episodes 
within a single very bright GRB, may help distinguish between the different 
possible causes for polarization and teach us about the dominant emission 
mechanism, the jet structure, or the magnetic field configuration within the 
GRB outfiow. 



11.10 Light curves for oflT-axis viewing angles 

The observed flux density = dE / dAdudt is the energy per unit area, 
frequency and time in the direction hd normal to a surface area dA (at the 
detector). The differential contribution to it is dF^{hd) = Ii,{n) cos 9 sdd^sd 
where cos Osd = n-fid, n being the local direction from the relevant part of the 
source to the observer (or detector), dQ,sd = dcpsdd cos Osd is the differential 
solid angle sustained by the contributing portion of the source as viewed 
from the observer, and luin) = dE / dAdQdudt is the speciflc intensity (the 
energy per unit area, time and frequency of radiation directed within a small 
solid angle d^l centered on the direction h). Since in practice almost always 
Osd ^ 1) one can simplify the expression for Fi, by approximating cos Ogd ~ 1- 
Furthermore, we have dQgd = dS±/d\ where dA{z) is the angular distance 
to the source, and dS± is the differential area in the plane of the sky (normal 

t A significant change in the brightness of the jet on an angular scale of ~ l/F around the line 
of sight can mimic a line of sight near the edge of a uniform jet, and produce a similarly high 
degree of polarization. 
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to n) sustained by the source, so that dFy = Ii,dS_\_/d\. Let us denote with a 
subscript z quantities measured at the cosmological frame of the source. For 
an optically thin source 1^^ = J ju^dsz, where j^,^ = dEz/dVzdi^zdi^zdtz is 
the emitted energy per unit volume, solid angle, frequency and time, while 
dsz is the differential path length along the trajectory of a photon that 
reaches the observer at the time tobs when is measured. Since Ivl^^-, 
ivjv^ and dsjv are Lorentz invariant (Rybicki & Lightman 1979), we have 
ly = {v/uzfly^ = (1 + S ju^dsz- In addition, since dA = di/ii + z)'^, 
where dL{z) is the luminosity distance, then dFi, = ju^dVz (1 + z)/d\, where 
dVz = dS±dsz is the volume element in the source cosmological frame. Here 
ji,^ = [r(l— n-/3)]^^ j'^i is measured in the source frame, while j'^, is measured 
in the (comoving) rest frame of the emitting material, which moves at a 
velocity /3c in the source frame. Altogether, this gives 

where u' = {1 + z)T{l — n ■ (3)v, and tz is the coordinate time at the source's 
cosmological frame. Since d'^x = dtzdVz = dtzdS±dsz = dtzdS^ds' {vz/^') = 
dtzdV'/T{l - n • /3) and AnjldV = dL'^, = ATT{dE' /dQ.'du'dt') is the differ- 
ential of the isotropic equivalent spectral luminosity in the comoving frame, 
eq. 111.81 can be rewritten as 

(1 + ^) f X *obs n-f\f dL'^, 



A7Tdl{z)J ' \' 1 + z c ; ; r3(l-n-/3)3 ■ 

(11.9) 

This result can be intuitively understood as follows. For a point source of 
luminosity L we have F = L/ATTd\ by the definition of the luminosity dis- 
tance, and thus F^, = dF/dv = {dL/duz){dvz/dv)/ATTd\ = L„^{1 + z)/A'Kd\, 
where L^^ = Ai^idEz/ dVlzduzdtz) = {dEz/ dE'){dVl' / dnz){dv' dt' / diyzdtz)L'^, 
= {^z/^')^L'^i = L'j^,/r3(l — n-/3)3 and we have used the relations dQ' /dQz = 
= iyzl^')'^ and dv'dt' = duzdtz derived in Sect. 111.21 This is the basic 
result for a point source, and the flux from a source of finite size can be ob- 
tained by dividing it into a large number of small regions that may be treated 
as individual point sources (in the sense that jy^ = [r(l - n ■ /3)]"^ j^, does 
not vary a lot within such a region) and summing over their contributions, 
as manifested in eq. 111.91 

In this section we are particularly interested in the dependence of the 
observed light curve, i^l, (tobs j "•) ) with viewing angle 0obs relative to the 
symmetry axis of the jet, that points at some direction nj, so that cos 0obs = 
h ■ fly For a jet that possesses such axial symmetry, the observed flux 
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depends only on ^obs and not on the azimuthal viewing angle (jjohs-, so that 
Fi, = -Fiy(tobs) ^obs)- In order to gain some intuition for how the afterglow 
light curve varies with viewing angle it is useful to first consider the simplest 
case of a point source moving at an angle of ^obs ^ 1 from the line of sight 
with r S> 1 where T oc R~"^^'^. In this case we can write 

where one uses the appropriate value of R{tohs, dohs) according to eq. 111.111 
For a given value of F, the Doppler factor P is roughly constant (~ F) 
for viewing angles within the "beaming cone" of half-opening angle 1/F 
around the direction of motion (i.e., for 0obs~l/r) and rapidly decreases 
as 2? oc 6~^^ for viewing angles outside the beaming cone (^obs > 1/r)- It 
is convenient to compare the observed flux density Fi, at different viewing 
angles ^obs and observed times tobs; that originate from the same emission 
radius R (Granot et al. 2002). In particular, Fy{Oohs > 0) with Fy{Oohs = 0), 

i^.(Us, ^obs) ^ a3F,/,(6tobs, 0) , a ^ ^^j^ « TTWUm^ ' 

^_ tobs(i?,0) 1 



" iobs(i?,^obs) " l + (l + rn)[F(i?)0obsF^''' ^^^'^^^ 

If the source decelerates, as expected during the afterglow (m > 0), then 
at early times F^obs ^ 1 and we have tobs ~ (1 + z)9'^^^R/2c oc R and 
a w (1 + m)6 [F(i?)0obs]~^ oc i?™ oc i^g. For an on-axis flux density of 

F,yiO,tohs) fx i^~'^t~^s, this implies that F^(6'obs, ^obs) oc _ 
If the jet is assumed to be non-expanding one can thus approximate F(i?) 
at the center of the jet by the expression for a spherical flow, for which 
m = 3 — k (see Sect. I11.5|) . so that for a stellar wind environment (A: = 2), 
the flux rises as t^^s'"^"^ ~ ^obs ~ *obs /3 ~ 1 and a ~ 1 — 1.5), while 
for a uniform external medium (A: = 0), the flux rises much more steeply 
as ^obs^'^"^" ~ *obs ~ ^obs- ^t late times when the beaming cone widens 
enough and engulfs the line of sight (F^obs < 1), the off-axis light curve 
approaches the on-axis one, since for F^obs ^ 1, a ~ 6 ~ 1 and, as a result, 

This simple model of a point source along the jet symmetry axis (called 
model 1 in Figure [Tl.9|) is instructive, as it captures much of the qualitative 
behavior of off-axis light curves. Nevertheless, it is important to also consider 
more realistic jet models. A useful model is a uniform jet of half-opening 
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Fig. 11.9. Afterglow light curves for different viewing angles (from Granot et al. 
2002): Left panel: comparing a simple analytic model featuring a point source 
along the jet axis (model 1; dashed lines) and a uniform jet with sharp edges that 
expands laterally at the sound speed in the jet comoving frame (model 2; solid 
lines). Right panel: comparing model 2 (inset) to the light curves calculated from 
a 2D hydrodynamical simulation of an initially uniform jet with sharp edges. 



angle 9j, with an initial value Oq, which expands laterally at some speed /SgC 
in its own, comoving rest frame (where /3s = corresponds to no lateral 
expansion and /3s ~ 1 corresponds to relativistic expansion in the comoving 
frame). For model 2 in Figure 111.91 /3sC is taken to be the local sound 
speed, which implies /3s ~ l/\/3 during the relativistic stage. The dynamics 
are given by conservation of energy and rest mass. The light curves for 
model 2 differ considerably from those for model 1 for ^o^s ~ at early 
times. This is because in model 2, for such viewing angles, the observed 
flux at early times is dominated by material along the line of sight, whose 
beaming cone encompasses the line of sight from the very beginning, while 
in model 1 there is no emission along the line of sight, and the emission 
from the point source along the jet axis is strongly beamed away from the 
observer. Moreover, in model 2 the light curves for ^obs ^ Eire very similar 
to those for ^obs = 0, since the observed emission is mainly from an angle 
of 1/r around the line of sight, and for such viewing angles this region is 
mostly (or even totally) occupied by the uniform jet, resulting in very small 
differences in the observed flux density. 

We note that model 1 can be made somewhat more realistic by placing the 
point source at the edge of the jet at the point P closest to the line of sight 
for viewing angles outside of the jet aperture, i.e., using 6 = max(0, 6'obs — ^j) 
instead of ^ = ^obs for the angle between the point source and the line of 
sight. In model 2, for viewing angles I/Fq <^ O^hs — ^ at early times 
when F(0obs — 6j) ^ I the emission is dominated by the part of the jet that 
is closest to the line of sight. Such emission arises from the area within the 
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jet where the Doppler factor is close to its highest value (at point P), which 
is over a solid angle of the order of ~ (^obs ~ ^j)'^- Therefore, if there were 
no lateral spreading, the initial rise in the flux would be similar to that for 
a point source at a fixed angle from the line of sight (which was discussed 
earlier). If, however, the jet spreads laterally, ^obs — d] would decrease in 
time, but the effect of the reduction in the fraction of the jet contributing 
to the observed emission is overwhelmed by the more rapid increase in flux 
due to the faster rate at which the beaming cone of the emitting material 
approaches the line of sight. 

For the prompt emission, the ratio of -Epeak ^ the photon energy at which 
the vFi, spectrum peaks - between lines of sight outside the jet (^obs > ^o) 
and those within the jet (which are similar to ^obs = 0) is a « + 
r^(0obs — ^o)^] ~ [r(^obs — ^o)]~^- The ratio of their apparent isotropic 
equivalent energy, E^^iso, which is the ratio of their observed fluence, is ~ 
for viewing angles that are relatively close to the edge of the jet, 1/Tq < 
^obs — ^0 < ^0 (Eichler &: Levinson 2004), since in this case only a fraction 
~ (^obs — ^o)^/^o ^/"^ of the jet (in the part closest to the line of sight) 
contributes significantly to the observed emission, resulting is a suppression 
by a factor of ~ 1/a relative to the result of for a point source. For 
a roughly standard -Epcak at viewing angles within the jet, this reproduces 

1/2 

the observed correlation -Epoak oc E^\^^, reported by Amati et al. (2002), 
Lloyd- Ronning & Ramirez-Ruiz (2002), and subsequently Lamb et al. (2005) 
using data from BeppoSAX, BATSE and HETE-2, respectively. For larger 
viewing angles, 6'obs^26'o, we have £'7,iso(6'obs)/-E'7,iso(0) ~ (r6'o)^[r(6'obs - 
9o)r^ ~ (^6lo)"^(6lobs/6'o)"^ which scales as ~ (r6lobs)"^- This is similar 
to the expected change of -E-y,iso with viewing angle for a point source, since 
for such viewing angles all of the jet contributes similarly to the observed 
emission, and it may be reasonably approximated as a point source. 

A more accurate description of the jet dynamics can be achieved with 
hydrodynamic simulations. The resulting light curves for an initially uniform 
jet with sharp edges, for different viewing angles, are shown in the right panel 
of Figure [TL9l (model 3). The initially sharp edges of the jet quickly become 
smoother and the jet becomes non-uniform, especially near the edges, were 
there is a sharp decrease in the Lorentz factor and energy per solid angle, 
and the velocity is not in the radial direction. This results in a much larger 
contribution to the observed flux at early times for large viewing angles 

t This empirical correlation has inspired many theoretical investigations arguing for its origin 
despite the many debates on its observational validity. The reader is referred to Gehrels et al. 
(2009) for a review. 



34 



J. Granot & E. Ramirez- Ruiz 



-3 



-3 



-3 



10' 



t (days) 



10' 



10 



12. 10" 



10"' 

10^ 
10" 
10" 
10"' 
10"' 
10"' 
10"' 

10^ 
10' 
10 



Hydrodynamic simulation of an 

initially uniform siiarp edged jet 




uniform jet with siiarp edges 




Gaussian jet: e, Tq-I °c e-^ /^e^ 




t (days) 



Fig. 11.10. Afterglow light curves for different jet structures, dynamics, and viewing 
angles (from Eichler & Granot 2006). The top panel is from an initially uniform jet 
with sharp edges whose evolution is calculated using a hydrodynamic simulation, 
the second panel is for a uniform jet with sharp edges, and the two bottom panels 
are for a Gaussian jet, in energy per solid angle, and either a Gaussian or a uniform 
initial Lorentz factor. The viewing angles are ^obs/^o = 0, 0.5, 1, 1.5, 2, 3, 4, 5, 
where is the (initial) half-opening angle for the uniform jet (two top panels) and 
the core angle [Oc) for the Gaussian jet (two bottom panels). 



outside the jet initial aperture, which in turn gives a slower rise in flux 
compared to a perfectly uniform jet with sharp edges (such as model 2). 

The light curves for different viewing angles depend not only on the jet 
dynamics but also on its initial angular profile, and in particular on Tq{6) 
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and S{9) (where £ is the jet energy per soHd angle). Figure [11.101 compares 
the hght curves for different viewing angles, between different jet structures. 
The first panel from the top is for model 3 described above, the second panel 
is for a uniform jet with no lateral expansion, while the last two panels are 
for a Gaussian £{0), with either a Gaussian or a constant Tq{9). For a 
Gaussian jet, if both £{6) and Tq{9) have a Gaussian profile (corresponding 
to a constant rest mass per solid angle in the outflow), then the afterglow 
light curves are rather similar to those for a uniform jet (Kumar & Granot 
2003). If, on the other hand, S{6) is Gaussian while ro(^) = const, then 
the light curves for off-axis viewing angles (i.e., outside the core of the jet) 
have a much higher flux at early times, compared to a Gaussian Tq{6) or a 
uniform jet (see the bottom two panels of Figure [Tl.lOp . due to a dominant 
contribution from the emitting material along the line of sight which has an 
early deceleration time in this case (Granot et al. 2005). 

11.11 Unification schemes: implications of viewing angle effects 

The appearance of GRBs depends so strongly on their orientation relative 
to us that our current classification schemes might easily be dominated by 
our random viewing angles rather than by more interesting (intrinsic) phys- 
ical properties. These inherently highly anisotropic GRB models imply a 
radically different GRB appearance at different viewing angles. In practice, 
GRBs of different orientations will thus likely be assigned to different classes. 
Unification of these fundamentally identical but apparently disparate classes 
is an essential stepping stone to studying the underlying physical proper- 
ties of GRBs. The ultimate goal is to discover which are the fundamentally 
important characteristics of GRBs - e.g., black hole mass, black hole spin, 
accretion rate - and how they govern the formation of jets, and the produc- 
tion of radiation. In this section we thus critically examine GRB off-axis 
models and contrast them with the afterglow observations of X-ray flashes 
and sub-luminous long GRBs. Since afterglow data in these cases are too 
sparse and insufficient to derive meaningful constraints on the overall pop- 
ulation, we concentrate our efforts on a few well-monitored examples. 

Empirical classification of GRBs: GRBs traditionally have been as- 
signed to different classes based on their duration. On the basis of this 
criterion, there are two classes of GRBs - short and long - dividing at ~ 2 s 
duration (Kouveliotou et al. 1993). GRBs have also been classified accord- 
ing to their spectral properties, albeit less successfully. In particular, bursts 
with lower photon energies -Epeak at which the vFy spectrum peaks have 
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been dubbed X-ray flashes (XRFs) based on observations by BeppoSAX, 
BATSE, and HETE-2 (Heise et al. 2001, Barraud et al. 2003, Kippen et al. 
2003, Lamb et al. 2004, Sakamoto et al. 2005). XRFs, with durations rang- 
ing from several seconds to a few minutes and a distribution on the sky 
consistent with being isotropic, are similar to long duration (>2s) GRBs. 
In addition to XRFs, the empirical classification of variable X-ray transients 
observed by HETE-2 had been expanded to include an intermediate class of 
events known as X-ray rich GRBs (XRGRBs). The spectrum of XRGRBs 
and XRFs is similar to that of long GRBs (Sakamoto et al. 2005) except 
for their lower values of -Epeak and of -E-y,iso (their X-ray and/or gamma-ray 
energy output assuming isotropic emission). In all respects XRFs, XRGRBs 
and long GRBs seem to form a smooth continuum. 

Many different models have been proposed for XRFs, most of which try 
to incorporate them into a unified scenario with GRBs. These models in- 
clude high-redshift GRBs (Heise et al. 2001), dirty (Iow-Fq) fireballs (Der- 
mer et al. 1999, Huang et al. 2002, Zhang et al. 2003), regular GRBs viewed 
off-axis (Yamazaki et al. 2002, Dado et al. 2004, Kouvehotou et al. 2004), 
photosphere-dominated emission (Drenkhahn 2002, Ramirez-Ruiz & Lloyd- 
Ronning 2002, Meszaros et al. 2002), weak internal shocks (low variability, 
AF < F; Zhang & Meszaros 2002, Barraud et al. 2003, Mochkovitch et al. 
2004), and large viewing angles in a structured (Lamb et al. 2005) or quasi- 
universal (Zhang et al. 2004) jet. Most of these models mainly aim at ex- 
plaining the low values of .Epeak in XRFs and do not address their expected 
afterglow properties. The afterglow evolution alone can, however, serve as 
a powerful test for XRF models. In fact, most of the models discussed for 
XRFs and XRGRBs have at least one major flaw in common: they do not 
naturally produce the very flat afterglow light curve seen at early times. In 
what follows we thus concentrate only on the class of models that naturally 
produce such light curves, that is, a roughly uniform jet with sufficiently 
sharp edges viewed from outside of the jet core. This class of models has 
been discussed quantitatively in Sect. 111.101 Since most afterglow data of 
XRFs and XRGRBs are too sparse, here we critically examine the role of 
off-axis models and contrast them with afterglow observations of XRGRB 
041006 and GRB 031203 for which the afterglow light curves are reason- 
ably well monitored from sufficiently early times. The reader is referred to 
Swift observations of XRF 080330 (Guidorzi et al. 2009) and GRB 081028 
(Margutti et al. 2009) for more recent examples of XRFs or GRBs, which 
bear similar characteristics to the ones discussed here. 



Off axis jet models of XRGRBs and XRFs: The case of 041006. 
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XRGRB 041006 was detected by HETE-2 (Galassi et al. 2004). It had a 
fluence of 5 x 10~^ erg cm~^ in the 2 — 30 keV range and 7 x 10"^ erg cm~^ 
in the 30 — 400 keV range, corresponding to /x/7 ~ —0.15 which classifies 
it as an XRGRB. It has a redshift oi z = 0.716 (Fugazza et al. 2004), 
which for a fluence of / w 1.2 x 10~^ergcm^^ in the 2 — 400 keV range 
gives i?-),,iso ~ 1-6 X 10^^ erg. It had an observed peak photon energy of 
-^pelk = 63t5 keV, corresponding to ^^^g!^^ = 1091^^ keV. Figure [TLTT] shows 
an off-axis model yielding an acceptable fit to the to the optical and X-ray 
afterglow observations of XRGRB 041006, which is also consistent with the 
upper limits at radio and sub-mm wavelengths (Granot et al. 2005). From 
this analysis one can conclude that a successful model for the afterglow of 
XRGRB 041006 is that of a collimated, misaligned jet interacting with a 
stellar wind external medium of mass density pext = Ar~'^, where r is the 
distance from the central source. The parameter values used in this fit are: 
E = 1.0x10^1 erg, = yl/(5 x lO" gr cm^^) = 0.03, 6*0 = 3°, Oohs = 1.156lo, 
p = 2.2, ee = 0.1, and €3 = 0.001. 

The optical light curve is very flat at early times (q ~ at t ^ 1 hr, where 
Fi, oc t~"v~^) and becomes steeper after a few hours (q ~ 1-2), which is 
a little steeper than the decay index in X rays at a similar time (q ~ 1 at 
t ~ Iday). Also, the ratio of the flux in optical and X rays at t ~ 1 day 
implies a spectral index of /3 ~ 0.7 — 0.75 assuming a single power law 
between them. This suggests that the cooling break frequency fc is above 
the optical after 1 day. Since one requires very extreme parameters to get 
fc to the X-ray range after 1 day (even getting fc to be above the optical 
after a day requires relatively low values of and of the external density), 
it is most likely that z^c is between the optical and X-ray at 1 day, which can 
also explain the steeper temporal decay index in the optical (by Aa = 0.25) 
for a stellar wind environment {k = 2). This favors a wind medium over a 
uniform density one, since otherwise the flux in the optical will decay more 
slowly than in the X-ray (also by Aa = 0.25), which is contrary to what 
is observed for XRGRB 041006. At t>5days there is a flattening in the 
optical light curve, which is probably due to an underlying SN component 
or host. This explains why the observed flux is higher than that predicted 
by our narrow relativistic jet model. 

For a GRB jet with well-defined edges, both the prompt gamma-ray flu- 
ence and the peak of the spectrum drop very sharply outside the opening 
of the jet, as (or depending on 6'obs) and a, respectively. Therefore, 
the low S^,iso of XRGRB 041006 combined with Ey^^i^o = E/{1 - cos^o) ~ 
E{2/0q) k, 7.3 X 10^^ erg (which serves as a proxy for the E'-y^iso value that 
would have been measured by an on-axis observer) implies (-Ek,iso/-£'7,iso)^^^ ~ 
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Fig. 11.11. A tentative fit to the optical R-band (upper curve) and X-ray (0.5 — 
6 keV, lower curve) data of XRGRB 041006 (see Granot et al. 2005 and references 
therein). The ROTSE-IIIa points are shown with asterisk symbols since they are 
unfiltered, but they can still be treated as R-band observations within the measure- 
ment errors. We also added two lines to the X-ray data which indicate the edges 
of the 1 fj confidence interval for the temporal decay index, a = 1 ± 0.1, and cover 
the duration of the Chandra observation. The inset shows the predicted spectral 
slope, —(3 = dlogF^/dlogh', in the optical (upper curve) and in the X-ray (lower 
curve), together with the values inferred from observations. 



6.8 and T ~ (-E'k,iso/-£'7,iso)"^''^(^obs — ^'o) ^ ~ 330. This implies a (cosmologi- 
cal) rest frame on-axis £'peak(6'obs < Oq) ~ (-E'k,iso/£^7,iso)^/^-E'peak ~ 740 keV, 
which falls closely within the observed ii'peak ~ -E'7,iso relation. 

Off-axis jet models for sub-luminous GRBs: the case of 031203. 

At a relatively small distance, with a redshift of z = 0.1055 (Prochaska 
et al. 2004), GRB 031203 was also atypical in its isotropic equivalent gamma- 
ray energy output, with E'-y^iso ~ 10^° ergs (Sazonov et al. 2004). In fact 
its E^^isQ was intermediate between GRB 980425 and more typical bright 
GRBs with E^;iso ~ 10^^ - lO^'^ergs (Frail et al. 2001). The gamma-ray 
light curve was smooth and similar to GRB 980425, consisting of a single 
peak lasting about 20 s and a peak photon energy E'peak > 190 keV. Soon 
afterwards, an optical counterpart was identified and follow-up observations 
by several telescopes revealed a supernova, SN 20031w, with a spectrum 
very similar to that of SN 1998bw (e.g., Malesani et al. 2004). Subsequent 
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X-ray observations of GRB 031203 with XMM and Chandra identified an 
X-ray source coincident with the optical transient. The flux decay rate and 
the isotropic luminosity of the X-ray afterglow also ranked the event as 
intermediate between GRB 980425 and classical GRBs (Kouveliotou et al. 
2004). A very faint counterpart was also detected at centimeter wavelengths 
where it displayed a peak luminosity more than two orders of magnitude 
fainter than typical radio afterglows (Frail et al. 2003), but again comparable 
to that of GRB 980425. 

GRB 031203, or at least its gamma-ray luminosity directed at us, was 
certainly very weak. A straightforward interpretation might be that the 
GRB was deficient in all its emissions in all directions. This idea is compat- 
ible with the afterglow light curve at radio frequencies. However, when one 
combines the fact that a 20 s long GRB was observed, as well as an X-ray 
and infrared afterglow, the situation is more constrained and in fact more 
consistent with a model in which GRB 031203 was a typical powerful jetted 
GRB viewed off-axis (Ramirez- Ruiz et al. 2005). 

The off-axis jet interpretation for GRB 031203 requires the viewing an- 
gle to have been ^obs ~ 26*0 (Figure [TTTT2]) . A misaligned jet with a typical 
energy expanding into a stellar wind with properties similar to those of Wolf- 
Rayet stars is thus consistent with the observations, especially with the slow 
initial decline rates seen in both the X-ray and radio afterglow. One question 
that naturally arises is whether the observed gamma-ray flux of GRB 031203 
can be explained within the framework of this model. The low E^^^sa of 
GRB 031203 implies Q = 3.8°(S^,iso/10^° erg)-V8(£'.^j3 x Qvgf/^ 
(rT/50)"^/'^, where Ejet is the kinetic energy of the jet, and T = ^obs/^o ~ 1 
and the fiducial values were chosen to match those of GRB 031203, which 
were either observed (-E^^iso ~ 10^*^ erg) or inferred from the fit to its after- 
glow (6*0 ~ 3° — 5°, £^jet ~ 3 X 10^° erg, T 1). A consistent solution for 
both afterglow and prompt emission can thus be found if F ~ 50 and T ~ 1, 
which imply more typical values of ~ 2 MeV (given the observed value 
-Epeak ~ 190 keV) and E^, 

iso ~ 10^"^ ergs when observed on-axis (consistent 
within the intrinsic spread of the -Epeak ~ E^^i^o relation). These results are 
applicable in the present context provided only that one further condition 
is satisfied, namely, that the (on-axis) jetted outflow be optically thin to 
high-energy photons. For a burst with -Bpeak ~ 2 MeV, F must exceed ~ 50. 
We thus conclude that the observations, especially the slow initial decline 
rates seen in the X-ray afterglow, are more consistent with an off-axis model 
in which GRB 031203 was a much more powerful GRB seen at an angle of 

f This follows from the scaling of -E^.iso with a, here assumed to be oc 
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Fig. 11.12. Afterglow emission from a sharp edged uniform jet in GRB 031203. 
Light curves are calculated for various viewing angles 6'obs for a GRB with the 
standard parameters E-^^t = 3 x 10^° erg, p = 2.4, = 0.15, es ~ 0.02, 0o = 5°, 
and = (M/lO-5M0yr-i)(w^/lO3kms-i)-i 0.1. The data for GRB 031203 
can be reasonably fit by different sets of model parameters (i.e., the parameters 
cannot be uniquely determined by the data). For example, a sharp-edged jet with 
^0 = 3.5° seen at ^obs ~ 2.256*0 gives also a reasonably good description of the 
observations provided that Ce ~ 0.1 and es = 0.04. 

about two times the opening angle of the central jet. 

Unified scliemes for GRBs: The empirical classification scheme by which 
an event is tagged as a GRB, sub-luminous GRB, XRGRB or XRF is rather 
arbitrary. Therefore, there could be some cases where a jet that is viewed 
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on-axis (^obs < ^o) will be classified as an XRGRB or XRF instead of as a 
GRB, or the opposite case in which a jet viewed ofF-axis (^obs > ^o) niight 
be classified as a GRB instead of as an XRGRB or an XRF. A more physi- 
cally motivated classification would be according to the ratio of the viewing 
angle 6'obs and the jet half-opening angle instead of relying purely on 
spectral characteristics as in the present empirical scheme. Such a classifi- 
cation would, however, be much harder to implement as it is not a trivial 
task to accurately determine the viewing angle. The strongest constraints 
could thus be obtained from afterglow light curves of XRFs an XRGRBs 
that are well monitored from early times and at various frequencies (rang- 
ing from radio to X-rays). Current examples include Swift observations of 
XRF 080330 (Guidorzi et al. 2009) and GRB 081028 (Margutti et al. 2009) 

Similarly, the large statistical sample of GRBs and XRFs with redshift 
measurements will allow a reconstruction of the intrinsic luminosity function 
of the prompt emission. If GRBs, XRGRBs and XRFs are only a manifes- 
tation of the viewing angle for a structured, universal jet (whose wings are 
producing the XRFs), then no break would be expected in the luminosity 
function. On the other hand, if GRBs are the results of viewing angles that 
intersect the jet (whether structured or not), while XRFs and XRGRBs are 
ofF-axis events, then one would naturally expect a break in the luminosity 
function. Guetta et al. (2005) found that a luminosity function with a break 
is favored in order for the predicted rate of local bursts to be consistent with 
the observed rate. This also prevents the existence of an exceedingly large 
number of GRB remnants in the local Universe. 

The relative fraction of XRFs and XRGRBs to GRBs is also expected 
to be different in the various models. If indeed an XRF corresponds to 
r(^obs — 0o)^di. few and (^obs — ^o) ^ ^O; then the solid angle from which an 
XRF is seen scales as /r or as for a constant F (at a constant distance 
to the source), while the solid angle from which a GRB is seen scales as 6q. 
Therefore, the ratio of solid angles for GRBs and XRFs scales as and 
more GRBs compared to XRFs would be seen for larger ^o- As the distance 
to the source increases, XRFs could be detected only out to a smaller off- 
axis viewing angle, while most GRBs would still be bright enough to be 
detected out to reasonably large redshifts. Therefore, the ratio of GRBs to 
XRFs should increase with redshift. Finally, if the true energy E in the jet 
is roughly constant, then the maximal redshift out to which a GRB could 
be detected would decrease with Oq since -E'-y,iso ex 9^"^. This would increase 
the statistical weight of narrow jets in an observed sample, as they could be 
seen out to a larger volume. 

Finally, we conclude with a few possible implications of the off-axis model 
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hypothesis for XRFs and XRGRBs. For sufficiently large viewing angles 

outside the edge of the jet, one might expect some decrease in the variability 
of the prompt emission. This is because the width of an individual spike 
in the light curve scales as At oc 1/a while the peak photon energy and 
fluence scale as a and — a^, respectively. Since the interval between 
neighboring spikes in the light curve is typically comparable to the width of 
an individual spike, At, then if At increases significantly for large viewing 
angles this would cause at least some overlap between different pulses which 
would smear out some of the variability. Thus one might expect XRFs to 
be somewhat less variable than GRBs, at least on average, where a lower 
variability might be expected for lower values of Ep^ak and -E-y^iso- This 
may lead to a simple physical interpretation of the observed variability- 
luminosity relation in the prompt gamma-ray/X-ray emission (Fenimore &; 
Ramirez-Ruiz 2000, Reichart et al. 2001). 

Another possible signature of the off-axis model for XRFs and XRGRBs is 
in the reverse shock emission. If the reverse shock is at least mildly relativis- 
tic, then the optical flash emission would be less beamed than the prompt 
X-ray or gamma-ray emission, due to the deceleration of the ejecta by the 
passage of the reverse shock. This might cause the optical flash to be sup- 
pressed by a smaller factor relative to the gamma-ray emission, compared to 
the corresponding on-axis fluxes. Thus XRFs or XRGRBs might still show 
reasonably bright optical emission from the reverse shock, which might in 
some cases be almost as bright as for classical GRBs. Finally, XRFs and 
XRGRBs might also show a larger degree of polarization compared to GRBs. 

Our understanding of GRBs has come a long way since their discovery over 
40 years ago, but these enigmatic sources continue to offer major puzzles and 
challenges. As we have described, our rationalization of the principal phys- 
ical considerations for the prompt and afterglow radiation emanating from 
these objects combines some generally accepted features with some more 
speculative ingredients. When confronted with observations, it seems to ac- 
commodate their gross features but fails to provide us with a fully predictive 
theory - but then again, no such theory exists as of yet. What is more valu- 
able, though considerably harder to achieve, is to refine models like the ones 
advocated here to the point of making clearer quantitative predictions, and 
to assemble, assess and interpret observations so as to constrain or refute 
these theories. What we can hope of our present understanding is that it 
will assist us in such an endeavor. 
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